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Abstract 

^ ! 

£Lt We propose the Proxemic Floor Field Model as an extension of the 

Floor Field Model, which is one of the successful models describing the 
pedestrian dynamics. Proxemic Floor Field is the Floor Field which cor- 
responds to the effect of repulsion force between others. By introducing 
the Proxemic Floor Field and threshold, we investigate the process that 
pedestrians enter a certain area. The results of simulations are evaluated 
* 53 ■ by simple approximate analyses and newly introduced indices. The dif- 

ference in pedestrian behavior due to the disposition of the entrance is 
<""] , also confirmed, namely, the entrance in the corner of the area leads to the 

£1, long entrance time because of the obstruction by pedestrians settling on 

the boundary cells. 

1 Introduction 

Study on pedestrian dynamics has attracted many physicists over the past few 
years [IJ [2J [3J. We can see many collective phenomena induced from the com- 
plexity of the system of pedestrian flow. In order to describe these phenomena, 
some models have been suggested and have achieved success in some situations. 
Helbing et al. designed the social force model [H0 E], which represents pedes- 
trian movement using differential equations. Interaction between pedestrians is 
demonstrated by introducing "social force" - which acts between others. Using 
this model, we can reproduce phenomena such as arching at an exit, lane for- 
mation in counter flow, and oscillations of flow directions at narrow doors where 
pedestrians are crowding in both sides. 
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Cellular Automaton (CA) based models are also often used for their simplic- 
ity. In CA based models such as the lattice gas model [7] and floor field model [9] , 
time and space are discretized and simple update rules are introduced. These 
characteristics of CA based models are suited for numerical simulations. 

The floor field model is one of the major CA based models that describe 
the pedestrian dynamics stochastically on latticed cells. In Ref.0, two types 
of Floor Field (FF) which control hopping probability are introduced. One is 
the Static Floor Field (SFF) which controls the "destination of pedestrians". 
The other is the Dynamic Floor Field (DFF) which describes the characteristic 
that pedestrians tend to follow their predecessor. In addition, other studies 
introduced the effects of wall jT2] , force [12] and collisions Q3] among pedestrians 
etc., and many valuable results have been obtained. 

Although CA based models have been studied extensively for describing the 
evacuation process [9J [TU1 HJ and counter flow [7] |H1 H], etc., few studies have 
been done for the proxemic problem. Proxemics is the term introduced by 
Hall Pj)| referring to the study of set measurable distances between people as 
they interact. Hall proposed four types of distance range: Intimate distance, 
Personal distance, Social distance, and Public distance. Hall also summarized 
the effects of proxemics as "Like gravity, the influence of two bodies on each other 
is inversely proportional not only to the square of their distance but possibly even 
the cube of the distance between them. " This repulsive force has been introduced 
into a few CA based mo dels [T5]. 

In this paper, the characteristics of the proxemic behavior are integrated 
into the floor field model. We just introduced the Proxemic Floor Field (PFF), 
which describes the proxemic behavior of pedestrians, and the threshold, which 
describes the "motivation to move" dependent on the crowdedness of the neigh- 
borhood of each pedestrian. By using this model, we investigate inflow process 
, in which the proxemic effect is predominant. Inflow process is the situation 
where pedestrians enter a certain area and acquire their own space; this can 
be observed on elevators, trains, attractions in theme parks, and many other 
places. This process is universally seen in the situations where pedestrians tem- 
porarily stay in a limited area. The study of this process can be applied to the 
design of buildings, infrastructures, and so on. In such situations, after entering 
the area, pedestrians look for the location where they can keep their own space 
and are not intruded upon by other pedestrians. This movement is driven by 
the repulsion force between others. After they successfully acquire the location, 
they stops walking. We also have to consider the attractiveness of locations. 
For example, the space in front of the exit in the train carriage is valuable for 
passengers in a hurry who do not want to be obstructed during boarding, pas- 
sengers who carry bulky items, and so on. On the other hand, the place may be 
uncomfortable for passengers not in a hurry because many passengers intrude 
their space every time the train arrives at the stations. 

Though there are some CA based models which consider the "force" of 
pedestrians[13l[T7l[T8l[T9], they mainly focus on the local effect related to phys- 
ical contact and thus they are not suited for simulating the non-local repulsion 
dominant in the inflow process. Moreover, in the social distance modetfTE\ deal- 
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ing with the long range effect of social force, proxemic force does not affect 
walking pedestrians, so this model cannot simulate inflow process either. Thus, 
the Proxemic Floor Field model is a new model capable of simulating inflow 
process. Furthermore, it should be noted that, although they seem to be just 
an opposite processes, the inflow dynamics are more complicated than the evac- 
uation ones because the destination is not clear in the inflow process; therefore 
we have to set the driving force on each pedestrian. This is why the inflow 
process has been less studied, compared to evacuation process. Of course, this 
model can also be used for other situations, which adds the effect of proxemic 
force to previous CA based models. 

We evaluated simulation results of the time required for all the pedestrians 
to enter the room by mean field analyses and confirmed the difference due to the 
disposition of the entrance. We focus on not only the inflow process, but also the 
steady state; the state in which every pedestrian has stopped. We evaluated the 
steady state by introducing a couple of new indices to represent the disposition 
of the pedestrians. 

2 Model 

In this section we summarize the update rules of Proxemic Floor Field Model. 
The space is discretized into cells which can either be empty or occupied by one 
pedestrian. Each pedestrian can move to one of the unoccupied next-neighbor 
cells (or stay on the present cell) at each discrete time step t — > t + 1 depending 
on Floor Fields and threshold as explained below. 

In the simulation of the inflow process, the Static Floor Field S represents 
the attractiveness / unattractiveness of each cell for pedestrians. We can control 
the pedestrian's behavior such as "stay on a comfortable place (local maximum 
/ minimum of S)" , "approach the destination (gradient of £)", or "never stay 
on an uncomfortable place (local minimum / maximum of 5*)", by setting the 
appropriate Static Floor Field. It should be noted that this Static Floor Field is 
often set as the value dependent on the distance from the exit in the context of 
the evacuation process. The Dynamic Floor Field D represents the characteristic 
that pedestrians tend to follow their predecessors. This tendency is observed in 
crowded or emergency situations. This is implemented by virtual pheromone, 
which has its own dynamics of diffusion and decay controlled by the parameters 
a and 8. Each pedestrian drops a boson on the present cell as he/she walks, and 
this boson decays or diffuses with probability a and S respectively in each time 
step. The Dynamic Floor Field D is the number of bosons remaining on the 
cell at each time step. These two types of Floor Fields have been historically 
used as fundamental and essential components of the Floor Field Model. 

As explained below, transition probabilities of Floor Field Model are modi- 
fied by following components newly introduced in this paper. 
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2.1 Proxemic Floor Field 



Proxemic Floor Field represents the effect of repulsion which acts between oth- 
ers according to their distance. As mentioned in Ref. [32], repulsion for small 
distance is much greater than that for large distance. In this paper we define 
the Proxemic Floor Field P^ generated by each pedestrian k as fD, and Pij 
is gained by taking sum of each P-j for all the pedestrians. Here, r*j is the 
distance between fc-th pedestrian and cell-(i,j). This distance is defined to be 
1 for the neighboring cells and its own cell, while other cells are measured in 
the L2 norm (see Figure [1} , for the naturalness of the model. By adopting this 
distance, Proxemic Floor Field generated by a pedestrian is uniformly assigned 
to neighboring cells, and this PFF does not influence its transition probability. 
(See figs. Hand El) 
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Figure 1: Distances from 
the 1st pedestrian in the 
center cell. Distances 
between two neighboring 
cells and within the cell 
itself are defined to be 1. 
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Figure 2: Proxemic Floor 
Field generated by one 
pedestrian. 
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Figure 3: Proxemic Floor 
Field generated by two 
pedestrians. 



2.2 Threshold 

Pedestrians stop walking after they acquire a suitable position in trains, ele- 
vators, etc. Here, we introduce a threshold which is related to the decision 
making about walk or stop: when the threshold is larger than the gradient of 
Floor Fields, pedestrians tend to stay on their present cells. Threshold 8 is 
given by ([3]) for each pedestrian. The subscript O corresponds to its present cell 
and Wij denotes total Floor Field , i.e., the weighted sum of Floor Fields. The 
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magnitude of Wij controls the pedestrian movement (explained by transition 
probabilities in the next subsection). The greater the Wij, the lower threshold. 
This corresponds to the fact that pedestrians in crowded area tend to move to 
more attractive places (low threshold), and that pedestrians in spacious areas 
tend to ignore the small gradient of total Floor Field (high threshold). 



6 = e max exp{-k T W } (3) 
Wij = —kpDij + ksSij + k P P tj (4) 

where, 

< S>7nax £ [0, oo ) :the coefficient which controls the magnitude of threshold. 

fcr £ (0, oo) :the sensitivity parameter which controls the steepness of threshold 
against Wij. 

ks G [0, oo) :the sensitivity parameter which weights the effect of Static Floor 
Field 

kjj G [0, oo) :the sensitivity parameter which weights the effect of Dynamic 
Floor Field 

kp G [0,oo) :the sensitivity parameter which weights the effect of Proxemic 
Floor Field, which is newly introduced in this paper explained be- 
low. 



2.3 Basic update rules 

Pedestrians follow the rules below in the Proxemic Floor Field Model. Note 
that each Sij is given and = for all cells at time step t = 

1. Dij is modified by the rules of diffusion and decay. 

2. Each pedestrian chooses one target cell based on the transition probabili- 
ties pij determined by ([5]) in the order numbered (sequential update) . 

{iVexp (—Wij — @) : transition to the neighboring cells,-, 
N exp (—Wo) ■ transition to the present cell 

Here, N is the normalization factor. After each pedestrian moves, proce- 
dures below are followed. 

• Whenever the pedestrian moves to its adjacent cell, D at its origin 
cell is added by one. 

• Pij is modified according to its transition. 

In this paper, we adopt the sequential update and reject other types of 
update methods (i.e. parallel and random update) for the following reasons. 
Parallel update is not adopted for the simplicity of the model since it will add 
the step of the conflict resolution to it. Since pedestrians essentially avoid 
others in the situation discussed in this paper, we can ignore the effect of the 
conflict, and by doing so, we concentrate on the dynamics which comes from the 
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introduction of PFF. When we deal with the situation in which conflicts occur 
(which is not discussed in this paper), the Proxemic Floor Field can be easily 
integrated into the model with the parallel update and the conflict resolution 
rule as shown in [TU]. Moreover, we also do not adopt the random update for 
the formation of inflow rule explained in the next subsection. Since the inflow 
speed is described as the inflow probability, a change of update order leads to 
the inconsistency of the inflow flux. 

The numbering of the updating order is assigned to each pedestrian from 
head to tail on the initial line as shown in Figure 4. This numbering corresponds 
to the fact that the each pedestrian mainly takes care of and is affected by other 
pedestrians who have already entered the area in this situation. 

The transition probabilities to the neighboring cell can be replaced by 



by modifying normalization factor N without a substantial change. This ex- 
pression is explained as the following: high probabilities are allocated to the 
neighboring cells with relatively low Wij , which lead to relatively low Wij — Wo 
, i.e., gradient of total Floor Field. This formation relates the gradient of to- 
tal Floor Field to the priority of directions. The term " —9 " reduces moving 
probabilities. As explained in the previous subsection, pedestrians sometimes 
stop walking, ignoring the small gradient of total Floor Field. This behavior is 
described by introducing this "braking" term; namely, the transition probability 
to the present cell (i.e., the probability with which particles stay on the present 
cell and do not move to neighboring cells) gets relatively high according to the 
magnitude of the threshold. Thus, transition probabilities are determined by 
taking Wij and into account. 

2.4 Inflow probability 

We introduce inflow probability which describes the inflow flux of pedestrians 
through an entrance - given by ([7]). In this paper, we discuss the entrance one 
cell wide. A^ X 4 denotes the number of pedestrians in a certain area showed in 
Figure HJ This expression corresponds to the fact that the inflow flux depends 
on the crowdedness of the area in front of the entrance. As as the area gets 
crowded, the speed of pedestrian inflow shrinks. On the other hand, inflow has 
the maximum speed (a = 1) if the area is free. p cr corresponds to the density 
where the inflow starts to decelerate. 



Pij = N'exp{-(W ij -W )-®} 



(6) 





P3x4 = 



12 



(8) 
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Figure 4: In each time step, one pedestrian enter the area from the entrance cell 
indicated by the cell colored dark gray with probability a. 3x4 area is indicated 
by the area colored gray (left). Figure on the right indicates the layout of the 
simulation and initial position of pedestrians. The number on each pedestrian 
represents the updating order. 

3 Simulation 

By using the PFF Model, we simulate the inflow process, and investigate the 
characteristics of the phenomena. In the following subsections, the number of 
pedestrians is fixed to 25 and we consider the square area of L x L. The Static 
Floor Field is simply set as 

„ f +oo , (i, j) : Entrance 

S « = { l eise ® 

to avoid the situation that a certain pedestrian stops on the entrance cell, hin- 
dering the entrance of the following pedestrian. This assumption is natural since 
the pedestrians - who know the entrance of the following pedestrians, do not 
stop walking just after their entrance in the real situations. On the other hand, 
here, the Dynamic Floor Field is ignorable, since pedestrians have no reason to 
follow others (c.f. evacuation, counter flow, etc.), or rather avoid others. Thus, 
we can set ko as 0. Since we are not concerned with irrational movement of 
pedestrians in this paper, we set enough large fcp and Q ma x- This decreases 
the probabilities of transition to undesirable cells. Therefore, the parameters 
remain in the model are Q m ax/kp and fcy. We rewrite O m ax/kp as O m ax in 
the following for the clarity of discussion. 

3.1 Entrance position and corner effect 

Let us investigate the effect of disposition of entrance to the time required. Here, 
the time required is defined to be the time step when the last pedestrian enters 
the room. We set the entrance of various positions named corner /m- cell (see 
Figure [5]). In order to focus on the effect of entrance disposition, we consider 
the condition of Q m ax — 0. The parameter p cr is set to be 0.2 and we have 
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acquired the simulation results as shown in Figure [5] It must be noted that 
only 8 cells are considered in calculating inflow probability of corner entrance. 

When the area is large enough, the time required monotonically decreases as 
the entrance gets close to the center. Because pedestrians tend to be attracted 
to boundaries for their repulsion effect, the pedestrians on the boundary hinder 
the inflow, when the entrance is near the boundary {corner). On the other hand, 
when the area is small, the time required has local minimum at 2cell. The long 
time required in the condition of Scell can be explained as the following: during 
the inflow process, the corner cell of the area is surely occupied, and next cell 
of the corner cell (corresponds to the entrance cell of lcell) is not occupied for 
the strong repulsion force of the pedestrian occupying the corner. Therefore, 
the cell 2 cells away from the corner tends to be occupied, and hinder the inflow 
when the entrance is set on 3cell. 

This result suggests to us the importance of entrance position especially in 
designing small areas for pedestrians. 

3.2 Relationship between threshold parameters and the 
time required 

The value of Q ma x and kx determine the shape of threshold. Large kr gives 
steep decreasing threshold against Wo, and <d max is the coefficient of threshold. 
Low kt gives gradual threshold and therefore we get relatively high values for 
large Wo ■ Ideally rational pedestrians can be considered to ignore the gradient 
of Wo less than or equal to threshold, thus we can describe the their sensitivity 
against proxemic pressure. In the case of inflow, sensitive pedestrians (with low 
threshold) leads less entrance time because they do not stay just in front of 
entrance for their repulsion force, and the density of 3 x 4 area tend to be kept 
low. By contrast, insensitive pedestrians (with high threshold) take more time 
than sensitive ones because they stay in front of the entrance until the area gets 
crowded. We simulate the situation of inflow with center entrance varying & m ax 
and hr from 0.1 to 2.0. The result of simulation is shown as Figure UJ, |U We 
can see the time required decrease almost monotonically as the Q m ax increases 
or k t decreases. 
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Figure 5: Variety of entrance position. The entrance cell m cells away from the 
corner is named "m cell" 
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Figure 6: The time required against entrance position for various L in the case 
of p cr = 0.2 and & max = 0. 

4 Analyses 

4.1 The time required and mean field analysis 

The time that it takes all the pedestrians to enter a certain area is the quan- 
tity of great concern in the context of engineering. Here, we evaluate the time 
required by mean field analysis. The inflow described in ([7]) decelerated under 
the condition of congestion. In other words, when the decelerate effect is dom- 
inant, the area is crowded to some extent and we can then consider that the 
deposition of pedestrians is averaged for their repulsion forces. On the other 
hand, when the area is not crowded, since pedestrians enter the area smoothly 
enough(a = 1), mean field analysis also gives good approximation. By consider- 
ing that the density of pedestrians is averaged and each pedestrian spreads into 
the area immediately, we can evaluate the time required where k is a variable 
corresponding to the numbering of pedestrians. 



TaV£ ~ ? a ^)~? min {l,T^} (10) 
P k = ^ (11) 

The results of simulation and approximate calculation are shown in Figure 
[9l Here, the entrance is set in the center of the boundary (see Figure and 
Qmax — for the simplicity of the argument. When the area is large enough, 
all the pedestrian enter the area in the probability a = 1. Therefore, the time 
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Figure 7: The time required against 
Qmax and kr in case of L = 7, central 
entrance, and p cr = 0.2 



Figure 8: The time required against 
©max and kr in case of L = 13, central 
entrance, and p cr = 0.2 



required converges to the number of pedestrians. By contrast, the small area 
leads to crowdedness and inflow is decelerated, thus the time required increases. 
In Figure [SI we can see that the time required for entrance increase rapidly in 
the cases of small area. 

4.2 Evaluation of steady state 

In the simulation condition we set, the steady state of the system corresponds 
to the situation where all the pedestrian have stopped (and no pedestrian will 
move once the system becomes steady), since we prohibited irrational movement 
by adopting large enough kp. To consider the actual problems such as the 
disposition or width of the doors of train carriages, the disposition of standing 
passengers is highly concerned. Basically, uniform dispositions of pedestrians 
are desirable from the perspective of effective use of space. To evaluate the 
desirability of disposition of pedestrians, we introduce two indices. One is spatial 
efficiency index E defined as the following, 

E = J2 W o (12) 

k 

Where Wq is the total Floor Field of the present cell of k-th pedestrian and 
the sum is taken over all the pedestrians. Since the values of total Floor Field 
corresponds to the uncomfortable quality of the cells, the sum of Wq can be 
viewed as "total proxemic stress of the system" . Namely, low E means the area 
is utilized efficiently. 

The other index introduced here is unevenness U defined as (|13|) 

U = ^ ~P (Imin) logP {Imin) (13) 

I,„i„ 

p(l mi n) = #{k\l<k<M,l^ in = l min }/M (14) 

Here, ^ iTl is the minimum distance between k-th pedestrian and others, and 
M is the number of pedestrians. This U can be viewed as the information 
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Figure 9: The time required against the size of the area in case of p cr = 0.2, 0.4, 
and Q max = 



entropy of minimum distances, and the value of U corresponds to the disorder 
of disposition of pedestrians. When the disposition is perfectly uniform (all the 
pedestrians give the same l m in), U = 0. On the other hand, when the disposition 
is perfectly disordered (no group of pedestrians give the same l m in),U = logM 

By introducing these E and U, the characteristics of the inflow can be inves- 
tigated and they can be used to categorize the inflow situations in more general 
cases. For example - as a more advanced application - the change of area shape 
or the introduction of obstacles can be considered. 

We can say that low E and low U give the desirable situations, and the 
opposite denotes the worst situations. Let us investigate the possible pairs of 
these values in the inflow situation. We plotted possible points of final states 
induced by many patterns of inflow generated by constant inflow probability 
a = 0.5 (Figure [10]) . Here, we fix kx — 1.0 and investigate them for various 
Qmax- E increases as Q ma x increases large, which accounts for the situation 
that pedestrians stay in front of the entrance, and the space is not utilized 
efficiently. We can also see the scatter of E is especially small when Q ma x is 
small. Turning our eyes to U, it should be noted that the scatter of U is large 
when Qmax is small and decreases as Q ma x becomes large. On the other hand, 
the average of U has local maximum around Q max — 1-0, and takes small value 
for large Q max . 
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Figure 10: Possible final states plotted on E-U plane for various <d m ax from 0.0 
to 2.0 in the case of p cr — 0.2, L = 15, a = 0.5, and kr = 1.0 

5 Conclusion 

For simulating the pedestrian dynamics in which proxemic effect is dominant, we 
presented a extended Floor Field model by introducing PFF (Proxemic Floor 
Field) and threshold. In this paper, we focused on the dynamics of "Inflow 
process" as the proxemic limit of pedestrian dynamics. By conducting some 
simulations and mean field analyses, we showed some meaningful results. In 
subsection 13.11 we not only found that the entrance on the corner cell causes 
long entrance time, but also that the local minimum of entrance time exists 
when the area is small. In subsection 13.21 the monotonous relationships be- 
tween the entrance time, and /ct, @max were confirmed, while in subsection 14. 11 
mean field analysis roughly reconstructed the monotonous relationship between 
the entrance time and the area size L. Finally, subsection 14.21 showed the char- 
acteristics of disposition of pedestrians at steady state by introducing newly 
introduced indices (i.e. energy E and unevenness U). Pedestrians with low 
threshold can lead to the final state of low E and scattered U. By contrast, 
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pedestrians with high threshold can lead to high E and more concentrated U. 

In this paper, we proposed the model which is capable of wide application 
and investigated the basic characteristics of it. The problems of wider entrance, 
appropriate setting of Sij conforming to the actual situations, and the effect of 
area shape etc. should be investigated in the near future. 
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